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Abstract Traditionally, the Lighthill-Whitham-Richards (LWR) models for homo- 
geneous and inhomogeneous roads have been analyzed in flux-density space with 
the fundamental diagram of the flux-density relation. In this paper, we present a 
I— I new framework for analyzing the LWR model, especially the Riemann problem at a 

p^ linear boundary in which the upstream and downstream links are homogeneous and 

"^ initially carry uniform traffic. We first review the definitions of local supply and de- 

(-H mand functions and then introduce the so-called supply-demand diagram, on which 

a traffic state can be represented by its supply and demand, rather than as density 
and flux as on a fundamental diagram. It is well-known that the solutions to the 
Riemann problem at each link are self-similar with a stationary state, and that the 
wave on the link is determined by the stationary state and the initial state. In our new 
framework, there can also exist an interior state next to the linear boundary on each 
>>• link, which takes infinitesimal space, and admissible conditions for the upstream and 

downstream stationary and interior states can be derived in supply-demand space. 

With an entropy condition consistent with a local supply-demand method in interior 

■^ states, we show that the stationary states exist and are unique within the solution 

framework. We also develop a graphical scheme for solving the Riemann problem, 
and the results are shown to be consistent with those in the literature. We further 



in 

o 

f^ discuss asymptotic stationary states on an inhomogeneous ring road with arbitrary 

^— I initial conditions and demonstrate the existence of interior states with a numerical 

►^ example. The framework developed in this study is simpler than existing ones and 

can be extended for analyzing the traffic dynamics in general road networks. 



1 Introduction 

Essential to effective and efficient transportation control, management, and plan- 
ning strategies is a better understanding of the evolution of traffic dynamics on a 
road network; i.e., the formation, propagation, and dissipation of traffic queues. The 
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seminal work by ( [LightMll and Whitham|[T955] l and (| Richards|[T956l l (LWR) at- 
tempts to study traffic dynamics with respect to aggregate values such as density p, 
speed V, and flux q. Based on a continuous version of traffic conservation 

and an assumption about the fundamental diagram of the flux-density relation q = 
Q{p), the LWR model of a homogeneous road link can be written as 

The corresponding speed-density relation is v = V{p) = Q{p)/p. Here the max- 
imum or jam density is denoted by Pj; i.e., p e [0,Pj]. Usually, V{p) is a non- 
increasing function of traffic density, v/ = V'(O) is the free flow speed, V(p/) =0, 
and q = Q{p) is unimodal with maximum flux or capacity C — Q{pc) where Pc is 
the critical density. Finally, traffic states with density higher than Pc are congested or 
over-critical, and those with density lower than Pc are free flowing or under-critical. 
Compared with microscopic traffic flow models (e.g. Gazis et aL]|1961 Nagel 



and Schreckenberg 1992) the LWR model can be used to analyze traffic evolution at 



the aggregate level with shock and rarefaction waves. With its analytical power and 
simplicity, the LWR theory has been extended for studying traffic dynamics in more 



general transportation networks. For examples, |Daganzo| ( |1997 ") proposed a traffic 
flow model for freeways with special lanes and high-occupancy vehicles with a two- 



regime fundamental diagram, and Wong and Wong ( 2002 1 proposed a multi-class 
model for heterogeneous drivers. 

In this paper, we are interested in the LWR model for a road with bottlenecks, 
where traffic characteristics such as free flow speed, jam density, the number of 
lanes, and capacity may be different for different locations. In other words, the 
fundamental diagram q — Q{x,p) depends on location. Such a road link is called 
inhomogeneous and the corresponding inhomogeneous LWR model can be written 
as 

^ + ^%^^0. (3) 

dt dx 

In order to understand the fundamental properties of equation ([3]), we usually ana- 
lyze its Riemann problem at x = 0. Hereafter, we will refer to the upstream branch 
as link 1, the downstream branch as link 2, and x = as a linear boundary. In the 
Riemann problem, links 1 and 2 are both homogeneous and initially carry uniform 
traffic. That is, 

(?(xp)-|2'(^)'-'<°' (4) 

^^'''^^"ie2(p),x>o, ^^> 

and 
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p(x,r = 0) = 



Pi,x<Q, 
P2,x>0. 



(5) 



Since ( |Mochon||1987| l, there have been many analytical and numerical studies 
related to the inhomogeneous LWR model in the literature. Roughly speaking, there 
have been two types of methods for solving the Riemann problem of inhomoge- 
neous LWR model n In the first type, the inhomogeneous LWR model can be ana- 
lyzed as a non-strictly hyperbolic conservation law (Isaacson and Temple 1992 Lin 



et al. II 1995) |Jin and Zhang 2003a ) or as a hyperbolic conservation law with a dis- 



continuous flux function dGimse and Risebro|19"90 |Gimse|1993[ Klingenberg and 



Risebro^l995^;^Dieh^l995 1996b; ■Pie hl and Wallin|1996||biehl|1996a[ [Zhang and 



Liu 2003 J Burger et al. 2005, 2008_), and various numerical methods can be used 
( |Ba'ie et al'.||2002[ |Zhan'g and' Liu||2'o05a|bl |Zhang et"aL]|2006t [Herty et al.pOOT] ) 



In the second type, the self-similar waves of the Riemann solutions are separated 
into links 1 and 2 by introducing a stationary state for each link, and the wave on 
each link is determined by a new Riemann problem of the corresponding homo- 
geneous LWR model (Seguin and Vovelle 2003 Garavello et al.||2007 i. Here the 
stationary states are subject to admissible conditions as weU as certain entropy con- 
ditions. This solution framework was first proposed for solving Riemann problems 
at general junctions with more than one upstream and downstream link ( Holden andj 
|Risebro_1995| |Coclite et all2005j . In (Seguin and Vovelle_2003j , the method was 
introduced for solving the inhomogeneous LWR model, and the stationary states are 
solved for a specific example. In ( Garavello et al.|2007[ l, a more general approach 
was proposed for solving the stationary states with a singular map method. How- 
ever, all these existing methods solve the Riemann problem in flux-density space: 
the first type of method is tedious due to the need to analyze kinematic waves on 
both links at the same time, and the second type of method fails to present the en- 
tropy condition in a physically meaningful way. In addition, all existing methods do 
not account for interior states in stationary shock waves ( van Leer||1984 Bultelle 



et al.[ 1998 ) and cannot be easily extended for studying traffic dynamics in a road 
network ( |jin||2003] l. Note that, in this paper, we do not intend to study numerical 
solution methods for solving the inhomogeneous LWR model. 

In this paper, we present a new framework for analyzing the inhomogeneous 



LWR model. We also adopt the method of wave separation by ( Holden and Risebro 
I1995\ but introduce a stationary state and an interior state for each branch. Here 
stationary states are the self-similar states at the boundary, and interior states do not 
take any space in the continuous solution and only show up in the numerical so- 
lutions as observed in ( jvan Leer|1984 Bultelle et al.|1998 ^. Rather than using the 
fundamental diagram, we introduce a so-called supply-demand diagram and discuss 
the problem in supply-demand space. After deriving admissible solutions for up- 
stream and downstream stationary and interior states in supply-demand space, we 



I 

In ( Daganzo 



2006 ; the inhomogeneous LWR model is solved in the space of cumulative number 



of vehicles as a calculus of variations problem, and the existence of its solution is proved for 
road links with point bottlenecks. However, the wave solutions of the Riemann problem are not 
explicitly discussed. 
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introduce an entropy condition based on the discrete supply-demand method ( |Da-| 



ganzo|[T995a Lebacquel|1996 i. We then prove that stationary states exist and are 
unique for given upstream demand and downstream supply, and interior states exist 
but may not be unique. Further we compare the Riemann solutions obtained by the 
new method with those in the literature for both the homogeneous and inhomoge- 
neous LWR models. We also apply the new framework for analyzing asymptotic 
stationary states on an inhomogeneous ring road and demonstrate the existence of 
interior states with numerical examples. 

The rest of the paper is organized as follows. In Section 2, we review the defini- 
tions of the supply and demand functions and the discrete supply-demand method 
for computing boundary fluxes. In Section 3, we introduce the supply-demand dia- 
grams and the structure of the solutions to the Riemann problem of the inhomoge- 
neous LWR model in supply-demand space. In section 4, we derive the admissible 
conditions for stationary and interior states in supply-demand space and an entropy 
condition consistent with the local supply-demand method in interior states. In Sec- 
tion 5, we solve the Riemann problem for both the homogeneous and inhomoge- 
neous LWR models and present a graphical solution scheme. In Section 6, we an- 
alyze asymptotic stationary states on an inhomogeneous ring road and demonstrate 
the existence of interior states with numerical solutions. In Section 7, we conclude 
our study with a discussion of future directions. 



2 Review of the supply-demand functions and methods 

2.1 Review of Engquist-Osher functions and the Godunov method 
for convex conservation laws 

For the original LWR model (|2|, assuming that k = Pc— p, we obtain a hyperbolic 
conservation law in ^ e [—p^Pj — Pc] as follows 

''•^m-o, (6) 



dt d 



X 



where /(A:) = C Q{Pc — k) is convex when Q{p) is concave, since J\ — 
-^Mfil > 0. Here /(O) = 0. Moreover, if ^ = Q{p) is the Greenshields funda- 



mental diagram (Greenshields|I935l, then (J6]) is Burgers' equation 



For the nonlinear equation (|6|l, we usually have to resort to numerical solutions 
for general initial and boundary conditions. After dividing the time duration into a 
number of time intervals of At and splitting the road link into a number of cells 
of width Ax, the finite difference equation in conservation form can be written as 
foUows ( [Colella and Puckett|2004[|LeVeque|2002l ): 
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r.^'+i 



At 
Ax 



( fi* _ fj* 

yi+\/2 Ji- 



1/2 



(7) 



where k\ is the average value of fc(x, f ) in cell / between x,- — ^ and jc,- + ^ at time 
step j, and //*, ,, is the flux through the boundary jci + 4^ between time steps j and 



'i+i/i 



j + 1. Here Ax and At have to satisfy the so-called CFL condition (Courant et al. 
[T928] l. 

For a hyperbolic conservation law ( 6|, the following functions were first intro- 
duced by ( [Engquist and Osher|1980a|bl|l98T] |Osher and Solomon] 1982| 

gW=/(max{^,0}) = |^|^)'i[J|U, 

k rk 

Xis)f'{s)ds^ / max{/'(.),0}fl'., 
Jo 

(l-ZW)/'W^-^= / min{f'{s),0}ds, 
Jo 



(8) 



(9) 



where j(A:) equals 1 iff /'(fc) > and equals otherwise. Note that f{k) = g{k) + 
h{k) and f{k) = msix{g{k),h{k)}. Therefore, we can rewrite (J6]) in the following 
form: 



kt+g{k)jc + h{k),,^0, 
kr + [max{g{k),h{k)}l^Q. 



(10) 
(11) 



Further, based on these definitions, the following E-O flux is introduced ( [En 
gquist and Osher|1980b[|Osher.l984) 



f!-i/2 = 8{kli)+h{ki)=f{kl,)+ / min{/'(.),0}fl'. 



fiki) 



mnx{f'{s),0}ds 






That is, (|7]) can be written as 



kt' =ki- 



At 
Ax 



{h{ki^,}-h{ki)+g{ki)-giki, 



(12) 



(13) 



This can be considered as upwind method for ( 10 1, since g{k) is non-decreasing and 
h{k) non-increasing. 

In ( |vanLeer|1984[ l, the Godunov flux (Godunov 1959 1 for Burgers' equation was 
written as 
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/:^j/2=maxte(fc^i),Mfc/)}- (14) 

In ( |Osher|1984| l, a new formulation of the Godunov flux ( |Godunov||1959| l was in- 
troduced as 

For convex f{k) — max{g(A;),/i(fe)}, since g{k) and h{k) are monotonically increas- 
ing and decreasing respectively, this is equivalent to 

,* _ , max{min^^^^^^^,gW,min^,^^^^^^,/2(^)}, ^_, < kj 



Ji-\I2 ^ 



max{max^j^^^^^^;g(^),max^;^^^^^j/2(^)}, k]_^ > kj 
= max{g{kl,)Mki)}. 



That is, ( 15 I is equivalent to ([14]). However, it has been shown that (15 1 can also be 



applied to non-convex f{k). 



2.2 Review of supply and demand functions and Godunov methods 
for the LWR model 

For the LWR model (J2]l, we define the following functions 



z,(p) = e(,^„(p,p,),.{«p).;;p|p. 



\is)Q'is)ds^ rmiix{Q'{s),0}ds (16) 

Jo 

= C+ r{l-xis))Q'{s)ds = C+ rmm{Q'{s),0}ds, (17) 

Jo Jo 

where x{p) equals 1 iff 2'(p) > and equals otherwise. It is straightforward to 
show that D{p) = C — g{k) and S{p) =C h{k). Therefore, the Godunov method 
for (|2]i is equivalent to 

where the boundary flux can be written as ( |vanLeer|1984| l 

<j/2=min{D(p,^i),5(p/)}, (19) 
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or as ( |OsherfT984l i 



7* 



1/2 



™Vi<p<p/^^^^' ^'^i<P'' 



max 



pi<P<pU 



Qip),Pi <PU 



(20) 



From (20 1, we can see that ( 19 1 is valid as long as the fundamental diagram Q{p) is 



unimodal and may not be concave. 

In the transportation literature, ( Bui et al.]|1992 1 first applied ( 1 8 1 and ( 20 1 for 
solving the LWR model. In ( Daganzo|1995a[ l, a new finite difference form was pro- 
posed for the LWR model with a triangular or trapezoidal fundamental diagram 



Q{p)^mm{vfP,Vc{Pj~p),Q,„ax} 



(21) 



where v^ = Vf-r^^ is the absolute value of the shock wave speed in congested 
traffic. In the so-called cell transmission model (CTM), the space-time domain was 



discretized with a CFL number of 1; i.e.. Ax = VfAt, and the boundary flux in (18 i 
was written as 



}* 
1- 



1/2 



min{D(p/li),5(p/)}, 



(22) 



where D{p^)At = mm{Q„uaAt,nj} and S{p^)At = nnn{QmaxAt,'^,{N,„ax - «/)} 
are defined as "the maximum flows that can be sent and received by cell / in the 
interval between time steps j and j + 1", nj — pjAx = pjvfAt is the number of 
vehicles in cell / at time step j, and Nmax = PjAx — PjVfAt is the maximum number 



of vehicles in cell /. Hence the physical meaning of ( 22 1 is that the boundary flux 



is the minimum of the upstream sending flux and the downstream receiving flux. It 
can be shown that D{p) = D{p) and S{p) = S{p). Thus, (22 1 is equivalent to ( 19 1 
for a CFL number of 1 and triangular or trapezoidal fundamental diagrams. 

Following ( jLebacque| 1 996) , we refer to D(p) in ( [T6| ) and S{p) in (17 1 as the de- 
mand and supply functions respectively and call \\9) the discrete supply-demand 
method for computing fluxes. The physical interpretations of demand and sup- 
ply functions and the supply-demand method have formed the basis for extend- 
ing the supply-demand method for computing fluxes through various network junc- 
tions ( Daganzol|1995a[ [Lebacque]|1996 ', 'Jin"2003). For inhomogeneous roads, the 
extension is straightforward as follows (Daganzo 1995a Lebacque 1996 ): the de- 
mand and supply functions in (fl6| and (fT7| are location-dependent, D{x,p{x,t)) = 
Q{x,m\n{p{x,t),Pc[x)}) and 5(x,p(x,f)) = g(jc,max{p(x,f),Pt(jc)}), and the flux 
is still computed by the supply-demand method ( 1 9 1, q'*^ j .j = min{D,_ i {Pj_ ^),Si (p/ ) }, 
where D,_ i (p ) is the demand function in cell i -T^ and St (p ) is the supply function 
in cell /. It has been shown that the flux by the extended supply-demand method is 
still the Godunov flux [hn and Zhang|2003a[ Zhang and Liu|2003 1. 

In this study, based on the Godunov finite difference equation in ( fTS] ! and the 
supply-demand method in ( [T9| for computing boundary fluxes, we attempt to con- 
struct the convergent solution of &\ with discontinuous flux functions (Hh and initial 
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conditions ([5]l. In ( |Daganzo| 1 995b| l, the convergence, truncation error, and capture 
of shock waves were directly derived from the corresponding finite difference equa- 
tion of the homogeneous LWR model. Here, we attempt to present a new framework 
and find the solutions to the Riemann problem of the inhomogeneous LWR model 
at a linear boundary. 



3 Supply-demand diagrams and the structure of Riemann 
solutions 

In the literature, the inhomogeneous LWR model has been analyzed in flux-density 
space. Since a traffic state in flux-density space {p,q) can also be represented in 
supply-demand space as U = {D,S), in this study we will analyze the inhomoge- 
neous LWR model in supply-demand space. Furthermore, we present the structure 
of solutions to the Riemann problem of the inhomogeneous LWR model in supply- 
demand space. 



3.1 Supply-demand diagrams 

Corresponding to the fundamental diagram in flux-density space, a supply-demand 
diagram can be introduced in supply-demand space. In Figure [T|b), we draw a 
supply-demand diagram for the two fundamental diagrams in Figure [it a). On the 
dashed branch of the supply-demand diagram, traffic is under-critical (UC) and 
U = {D,C) with D < C; on the solid branch, traffic is over-critical (OC) and 
U = {C,S) with S <C. Compared with the fundamental diagram of a road section, 
the supply-demand diagram only considers capacity C and congestion level of traf- 
fic flow, but not other detailed characteristics such as critical density, jam density, or 
relationship between density and flux. That is, different fundamental diagrams can 
have the same supply-demand diagram, as long as they have the same capacity and 
are unimodal, and their critical densities, jam densities, or shapes are not relevant. 
However, there is always a one-to-one mapping between a given supply-demand di- 
agram and its corresponding fundamental diagram. That is, there exists a one-to-one 
mapping between (p.q) and {D,S). 



For the demand and supply functions in (16 1 and (17 1, we can see that D is non- 
decreasing with p and S non-increasing. Thus D <C, and S <C, and max{D, S} — 
C. In addition, D = 5 = C iff traffic is critical; D < S — C iff traffic is strictly under- 
critical (SUC); S < D = C iff traffic is strictly over-critical (SOC). Therefore, the 
state U = {D,S) is under-critical (UC), iff S ^C, or equivalently D < S. The state 
U = {D,S) is over-critical (OC), iff D = C, or equivalently S<D. 

From a state on the supply-demand diagram, we can obtain the correspond- 
ing flux q{U) = min{D,S} and capacity C = max{D,S}. However, we cannot tell 
the density from the supply-demand diagram, and the fundamental diagram is still 
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(b) 



Fig. 1 The fundamental diagram and its supply-demand diagram 



needed to compute the density p from {D, S). That is, p can be written as a function 

of (D, 5) as 



PiD,S) - . ^_i 



D-\D),D<S, 
S-\S), S<D, 



since D{p) and S{p) are invertible when the traffic is UC and OC respectively. Note 
that p is not a function of q. If we introduce the supply-demand ratio 7 = D/S, then 
q{D,S) = OT/n{7, 1/7} • C, and 

D-i(C7), 7<1, 



p(D,5)^/.(7).|,-.(^,;;';;i; 

where R{y) is an increasing function in 7 e [0,°°]. Here R{0) = 0, R{1) = Pc, 
and R{°a) = Pj. In this sense, p = R{y) can be considered as the inverse flux- 
density relationship. Similarly, v = V{p) ~ V{R{y)) is a non-increasing function 
in 7, y (0) = v/, and y(oo) = 0. 



3.2 The structure of solutions to the Riemann problem 



In supply-demand space, initial conditions (J5]l are equivalent to 

^ ' ^ \U2^{D2,S2),X>0. 



(23) 



The Riemann problem at a linear boundary is then equivalent to the Riemann prob- 
lem for (|3]l with initial conditions ( |23] l. 

Unlike existing studies of hyperbolic conservation laws with discontinuous flux 
functions, in which solutions to the Riemann problem have been constructed in p — ^ 



space, within the framework of wave separation developed in (Holden and Risebro 
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|1995| l, we construct the solutions to the Riemann problem for ^ with initial condi- 



tions (23 1 in supply-demand space. 



In solutions to the Riemann problem for ([3]) with initial conditions ( 23 1, a shock 



wave or a rarefaction wave could initiate on a link from the linear boundary x = 0, 
and traffic states on both links become asymptotically stationary after a long time. 
We denote the stationary state on the upstream link 1 by U^ and the stationary state 
on the downstream link 2 by f/j^. At the boundary, there can also exist interior states 
(van Leer| 1984 Bultelle et al.|1998)l, which take infinitesimal space. We denote the 



interior states on links 1 and 2 by f/j (0 , t) and U2 {0 ,t) respectively. The structure 
of Riemann solutions on upstream and downstream links are shown in Figure [2] 



f/^(o+.t) 

A— 



(a) (b) 

Fig. 2 Structure of Riemann solutions: (a) Upstream link 1; (b) Downstream link 2 

Then the kinematic wave on the upstream link 1 is the solution of the correspond- 
ing homogeneous LWR model with initial left and right conditions of Ui and U[ , 
respectively. Similarly, the kinematic wave on the downstream link 2 is the solution 
of the corresponding LWR model with initial left and right conditions of U2 and 
U2, respectively. Since the stationary and interior states are constant for f > 0, states 



on both links 1 and 2 are self-similar (Smoller||1983). That is, if stationary states 



exist and are unique, we have unique self-similar solutions for the Riemann prob- 
lem of ([3]l. In the following sections, we first derive necessary conditions for both 
stationary and interior states and then solve them in supply-demand space. 



4 Necessary conditions for the existence of stationary and 
interior states 

We denote ^1^2 (0,f) as the flux from link 1 to link 2. We first observe that the 
fluxes are determined by the stationary states: the asymptotic out- flux of link 1 is 
qi{Q~,t) = q{Uy), and the asymptotic in-flux of link 2 is ^'2(0+,?) = ^'(t/j^). Fur- 
thermore, from the conservation of traffic at the linear boundary, we have 

^1^2(0,0 = qi{Q-,t)^q2{0+,t)^q{U^)^q{U+). (24) 
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4.1 The admissible conditions for stationary states 



11 



As observed in (Holden and Risebro 1995 Coclite et al.||2005 i, the speed of the 
kinematic wave on an upstream link cannot be positive, and that on a downstream 
link cannot be negative. We have the following theorem. 



[Ci.Di) D, 




D 




Fig. 3 Admissible stationary states for the upstream link 1 



Theorem 1 (Admissible stationary states). For initial conditions in ( [25] ), station- 
ary states are admissible if and only if 



f/f = (Di,Ci)or(Ci,5^), 



(25) 



where 5j < D\, and 



U+^{C2,S2)or{D+,C2), (26) 

where D^ < 82- The results are illustrated in Figuresplandm 

Proof. When the initial state on link 1 is strictly under-critical; i.e., when Di < ^i = 
Ci, the admissible stationary state f/f is the same asUi = {Di,Ci) or strictly over- 
critical with U[^ = {Ci,S^), where S^ < D\. In this case, the Riemann problem 
for the LWR model on the upstream link 1 with upstream and downstream initial 
states U\ = {Di,Ci) and U[' has the following possible solutions: there is no wave 
when f/f = Ui; there is a backward traveling shock wave when U^ = {Ci,S^). In 
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(a) 



Ij (S2,C2) 



I If/2 = U} 



(c) 




Fig. 4 Admissible stationary states for the downstream link 2 



addition, we can verify that any stationary states not equal to L'^i and with S^ > Di 
will lead to forward traveling shock waves or rarefaction waves. Note that when 
U^ = (CijDi), the Riemann problem is solved by a zero shock wave, but f/^^ cannot 
be the stationary state by definition. When the initial state on link 1 is over-critical; 
i.e., when Si < D\ = C\, the admissible stationary state U^ is over-critical with 
f/j~ = (Ci,5]"), where 5^ <C\. In this case, the Riemann problem for the LWR 
model on the upstream link 1 with upstream and downstream initial states U\ = 
(Ci,5i) and U^ — {C\,S^) has the following possible solutions: there is no wave 
when S^ — Sy; there is a backward traveling shock wave when ^i > 5j" ; and there 
is a backward traveling rarefaction wave when S\ < S^ . Therefore, the stationary 
state is indeed admissible. In addition, we can verify that any strictly under-critical 
stationary states U^ will lead to forward traveling rarefaction waves and are not 
admissible. 

When the initial state on link 2 is under-critical; i.e., when D2 < S2 = C2, the 
admissible stationary state U2 is under-critical with U2 = {D2 ,€2), where D2 < 
C2. In this case, the Riemann problem for the LWR model on the downstream link 2 
with upstream and downstream initial states U2 = {Dj ,€2) and U2 = {D2,C2) has 
the following possible solutions: there is no wave when D2 = D2', there is a forward 
traveling shock wave when D2 < D2', there is a forward traveling rarefaction wave 
when D2 > D2- Therefore, the stationary state is indeed admissible. In addition, we 
can verify that any strictly over-critical stationary states U2 will lead to backward 
traveling rarefaction waves and are not admissible. When the initial state on link 2 is 
strictly over-critical; i.e., when S2 < D2 — C2, the admissible stationary state U2 is 
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the same as U2 — (C2,52) or strictly under-critical with U2 — {D2 ,€2), where D^ < 
^2. In this case, the Riemann problem for the LWR model on the downstream link 2 
with upstream and downstream initial states U2 andf/2 = {D2,C2) has the following 
possible solutions: there is no wave when U2 — U2', there is a forward traveling 
shock wave when U2 — {D2,C2)- In addition, we can verify that any stationary 
states not equal to U2 and with Dj" > S2 will lead to backward traveling shock 
waves or rarefaction waves. Note that when U2 = (52, C2), the Riemann problem is 
solved by a zero shock wave, but U2 cannot be the stationary state by definition. ■ 

Remark 1. Note that f/f = U\ and U2 = U2 are always admissible. In this case, 
the stationary states are the same as the corresponding initial states, and there are no 
waves. 

Remark 2. From the proof of Theorem [T] we can see that the types and traveling 
directions of waves on a homogeneous road can be solely determined by upstream 
demand and downstream supply and are not related to the shape of fundamental 
diagrams. This is why we are able to discuss the Riemann problem of the inhomo- 
geneous LWR model in supply-demand space. Note that, however, the wave speeds 
are related to the details of the flux-density relation Q{p)- 

Remarks. Then the out-flux ^i(0^,f) —q{U^) —mm{D^,S^} <Di and the in- 
flux ^2(0^,0 — q{U2) — minjDj ,5^} < ^2. That is, D\ is the maximum sending 
flux and ^2 is the maximum receiving flux in the sense of ( Daganzo| 1 994| [r995a I. 



Furthermore, ^i(0",f) = Di, iff U^ = (£>i,Ci), and iff f/f is UC; q\{0^ ,t) < Di 
iff t/f = (Ci,57) wifli 5f < Di, and iff f/f is SOC. Similarly, ^2(0+,?) = ^2, iff 
U^ = (C2,52), and iff U^ is OC; ?2(0+,0 < ^2 iff U^ = (£>2",C2) with dJ < ^2, 



and iff U^ is SUC. 



4.2 The admissible conditions for interior states 

The Riemann problem on link 1 with left and right initial conditions of U^^ and 
t/i (0^ , f ) cannot have negative waves. Otherwise, Ui (0^ , f ) will propagate upstream 
and violates the condition that it only exists at the boundary, but not anywhere up- 
stream. Similarly, the Riemann problem on link 2 with left and right initial con- 
ditions of U2{0^,t) and U2 cannot have positive waves. Therefore, interior states 
f/i(0^,f) and {/2(0+,f) should satisfy the following admissible conditions. 

Theorem 2 (Admissible interior states). For asymptotic stationary states Uy and 
U2 , interior states Ui{Q^ ,t) and U2{0^ ,t) are admissible if and only if 

^/, rn- / ^ - / ^*^i ' ^1^ ) = ^1 ' ^''''" s^ <D- =Ci 

uiyy ^') - \(^Di{0-,t),Si{0-,t)),whenDy <Si =Ci ^ ' 

where Si (0^ , — ^F' '^'^'^ 

lh(Q+ t\-l (^2+'^2) = U+, when D+ < S+ = C2 

t^2iu ,t) - I (£,2(0+,f),52(0+,f)), whenS+<D+^C2 ^^^ 
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it/i(0-,t) 

(Ci,flr) 



(a) 



(«f =fi(0-,t) 



(c) 




Fig. 5 Admissible interior states for upstream link 1 



where D2(0^,f) > 82- The results are illustrated in Figures\5\and\6\ 

Proof. The results can be verified with the observation that the Riemann solutions 
for the homogeneous LWR model of the upstream links cannot have negative waves, 
and those of the downstream links cannot have positive waves. ■ 

Remark 1. Note that U\{Q^ ,t) = Uy and C/2(0^,f) — U2 are always admissible. 
In this case, the interior states are the same as stationary states, and it is equivalent 
to saying that there are no interior states. 



4.3 An entropy condition for the local supply-demand method 



In addition to traffic conservation and admissible conditions, we introduce an en- 
tropy condition such that the boundary flux is always consistent with that by the 



supply-demand method (19 1 for the local interior states. At the boundary at x ~ 0, 
the immediate upstream state is Ui{Q^,t), and the immediate downstream state 
1/2(0^, f)- That is, the interior states have to satisfy the following entropy condi- 
tion 



^1^2(0,0 - mm{Di{Q-,t),S2{0+,t)}. 



(29) 



Note that the entropy condition ( [29] l is also equivalent to the following localized 
optimization problem 
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(a) 



, >f./+ 



D 



(c) 




Fig. 6 Admissible interior states for downstream link 2 



max{^i^2(0,0=?i(0 ,0=?2(0+,0} 
subject to 

qi{0-,t)<Di{0-,t), 
qiiO+j) < 52(0+, f). 

That is, the stationary and interior states are solutions of the optimization prob- 



lem in the domain defined by the traffic conservation condition (24 1, the admissible 



conditions for stationary states ( 25p6 i, and the admissible conditions for interior 
states ( 27p8 i. Optimization formulations of entropy conditions were also adopted 
in ( Ho den and Risebro|1995 Coclite et al.|2005 l, but in terms of stationary states 
in flux-density space rather than in supply-demand space as we have done here. 

In all of the necessary conditions above, we can see that the stationary and inte- 
rior states are independent of the upstream supply, ^i , and the downstream demand, 
D2- That is, the same upstream demand and downstream supply will yield the same 
solutions of stationary and interior states: when the upstream traffic is congested, its 
congestion level is not relevant to the stationary and interior states or the boundary 
flux; when the downstream traffic is free flow, its density is not relevant to the sta- 
tionary and interior states or the boundary flux. Note that, however, the types and 
speeds of waves on both links can be related to Si as shown in Figure [3|d) and D2 
as shown in Figure |4|d). 
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5 Solutions to the Riemann problem for the inhomogeneous 
LWR model 



For the Riemann problem of ([3]) with initial conditions ( 23 i, we first solve for the 



stationary and interior states that satisfy the traffic conservation condition ( 24 1, the 



admissible conditions for stationary states ( [25p6 l, the admissible conditions for 



interior states ( [27p8] l, and the entropy condition (|29|l. Then we compare them with 
existing solutions for both the homogeneous and inhomogeneous LWR models. 

Lemma 1. In the Riemann solutions, the boundary flux satisfies 

qi^2[0.t)^q{U^)=q{U+)=mm{D,,S2)- (30) 

Proof. From Theorem[T]and traffic conservation ( [24| , we have that 

We demonstrate that it is not possible for ^1^2(0,?) < min{Di,52}. Otherwise, we 
assume that ^1^2 (0,0 =qo <min{Di,52} < minjCi , Ci} - Since q{U^) =qo<Di, 



from (25 I we have that U^ = {Ci,qo). Further from (27i we have that L'^i(0^,f 
f/j" — (Ci,qo). Hence Di(0",f) =Ci. Similarly, since^f/j^) — qo < S2, from (26i 
we have that U2 = (^cQ)- Further fro m ([28| ) we have that f/2(0^,f) = UJ = 
{qo,C2)- Hence 52(0+;f) = €2- Then from (|29[) we have ^'1^2(0,?) = min{Ci;C2}, 
which contradicts the assumption that g'i^2(0,f) < min{Ci,C2}. Therefore we have 

m 



That is, the local optimal solution in ( 29 1 leads to a global optimal flux at a linear 



boundary, which satisfies the following optimization problem 

max{^i^2(0,f)} 
subject to 

^1^2(0,0 <Di, 
qi^2{0,t) < S2. 

Therefore, in the Godunov finite difference equation ( fTS) , the boundary flux at 
the first time step by (fT9]l is the same as the asymptotic flux, regardless of the time 



step size. That is, the discrete flux in ( 19 1 is consistent with the continuous flux in 



Theorem 3. The stationary states and interior states of the Riemann problem for 
pi) with initial conditions ( |23| ) are the following: 

1. When D\ < S2, we have unique stationary and interior states: U^ = U\ (0^ , t) = 
(£)i,Ci) andU} = f/2(0+,f) = (Z)i,C2); 

2. When D\ > S2, we have unique stationary and interior states: U^ = U\ (0^ , t) = 
(€1,82) andU+ = U2i0+,t) = (C2,52); 
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3. When D\ — S2, we have unique stationary states: t/j — (Di,Ci), and U2 — 
(C2,52); but interior states may not be unique: f/i(0^,f) = (Di,Ci) and D2{0^ ,t) > 
S2and 82(0+, t)>S2,orU2iO+,t) = {€2,82) and Si{0-,t)>Di and Di{0-,t)> 
Di. 



Proof. When Di < S2 < C2, ^ leads to ^1^2(0, = £>i. For link 1, since 
q{U^) — D\, from (25i we have that U^ — {Di,C\). Further from (27i we have 
thatf/i(0",f) = (Di(0"^,f),5i(0",f)) with5i(0",f) >£>i" =£>i. For link 2, since 
gr(f/+) = Di < 52, from (|26|l we have that U:^ = {Di,C2). Further from (|28i 



we have that 1/2(0+,?) = U^ = (Z)i,C2). Then from (|29J) we have ^1^2(0,?) = 
min{Di(0",f),C2} =£>!. Therefore, Di(0",f) =Di and5i(0",f) =Ci. 

When 52 < Di < Ci, (|30]) leads to q\^2{0,t) = 82- For link 1, since q{U^) = 

(Ci,52). Further from \ 2l\ 
(Ci,5'2). For link 2, since ^(C/2^) = ^2, from 326] 



f/i(o-,f) = t/r 



^2 < D\, from (25 I we have that U 



we have that 
we have that 
[/2+ = (C2,52). Further from (|28jl we have that t/2(0+,f) = (D2'(0^,f),52(0+,f)) 
with 1)2(0+, f) > 52. Then from"p9]l we have ^1^2(0,?) = min{Ci,52(0+,f)} = ^2. 
Therefore, 52(0+,f)= ^2 and Z)2(0+,f) = C2. 

WhenDi =^2, (|30| leads to g'i^2(0,0 =£)i ^ 52- For link 1, since q{U^) = D\, 



from ( 25 I we have that f/[ = (D 1 , Ci ) . Further from ( 27 1 we have that t/i (0 , f ) = 
(Di(0 



,f),5i(0",r)) with 5i(0",f) > D\ = D^. For link 2, since q{U:{) = S2, 
from (|26|l we have that t/j^ = (€2,52). Further from (|28[) we have that {/2(0+,f) 



(D2(0^),52(0+,f)) with Z)2(0+,f) > 52- Then fronr|29| we have ^1^2(0,?) = 
min{Di(0-,f),52(0+,f)}=Di =52- If Di(0-,f) = Dj == ^2, then 5i (0^,0 =Ci, 
and 52(0+, f) > 52- In this case, the interior state t/2(0+,f) may not be unique with 
■D2(0+,f) > 52 and 52(0+, f) > 52- Note that, when 52 = C2, the interior state is 



= 52, thenD2(0+,f) == C2, andDi(0",f) >Di. 
, f) may not be unique with 5i (0^,?) > D\ and 



£/2(0+,f)=f/2+-If52(0+,f)=£>l^ 

In this case, the interior state U\ (0" 
^i(0",0 > D\. Note that, when Di = C\, the interior state is U\{Q' ,t) = t/f . If 
both the upstream link 1 and the downstream link 2 have the same fundamental di- 
agram, this case corresponds to a stationary shock, and the interior state is the same 
as that in ( |vanLeer|1984^ . ■ 

Remark 1. From the theorem we can see that the stationary states always exist 
and are unique for the same pair of Di and 52- Thus, with given U\ and U2, we can 
find uruque kinematic waves on both links 1 and 2. Therefore, in the new solution 
framework, the solutions for the Riemann problem of the inhomogeneous LWR 
model always exist and are unique, although we may have multiple interior states. 



Remark 2. If the entropy condition (29i is replaced by (30i, we still have the 
same solutions t/j" and U2 ■ That is, if we do not consider possible interior states as 



in (Seguin and Vovelle|2003l|Garavello et al.|2007), then traffic conservation (|24l), 
I — 1 [i _ — ' V— K 

admissible conditions for stationary states ( 25p6 i, and the entropy condition ([30| 
will yield the same stationary state solutions. However, this simplified approach 
- which is what currently exists in the literature - does not yield the existence or 
properties of the interior states. 

Remark 3. When Di = S2, we have the following interior states that are different 
from the stationary states at jc = 0^ or ;c = 0+ . The interior state at jc = 0^ has to 
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satisfy ^(t/i(0 ,t)) >Di, and the interior state atx = 0+ has to satisfy q'(f/2(0^,f) > 

Si- 

In addition, we have the following conclusions concerning possible stationary 
states at a linear boundary. 

Corollary 1. When both links 1 and 2 reach asymptotic stationary states, they share 
the same flux q, and possible stationary states are the following: both links are UC 
with link 1 at {q,Ci) and link 2 at {q^Ci) where q — D\ < S2; both links are OC 
with link 1 at {C\,q) and link 2 at {C2,q) where q — S2 < Di; link 1 is UC at {q,C\) 
and link 2 OC at {C2,q) where q = D\ =82-11 is not possible that link 1 is SOC and 
link 2 sue 

Remark. The stationary states are stable in the sense that, when they are given as 
initial states, we obtain the same stationary states following TheoremIS] 



5.1 The homogeneous LWR model 



(a) 



(d) 



(D,.S,) = f',- = !/^ 



(/j = ur = Ut 



(b) 




r D 



(e) 



(f) 



Fig. 7 The Riemann problem for the LWR model: Stationary states in supply-demand diagrams 



For the original LWR model (|2]), the upstream link 1 and the downstream link 



2 have the same fundamental diagram. Therefore we have Ci = C2- In (Lebacque 



1996[ l, there are four scenarios for solutions to the Riemann problem. Here we re- 
organize them into the following six cases of initial conditions: 
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Case 1 Link 1 is SUC, and link 2 is UC. That is, Di < Si = C\ and D2 < S2 = €2- 
In this case, Di < 82- From Theorem pi we have that U[' — Ui{0^ ,t) — U2 — 
t^2(0^,f) = Ui and ^1^2(0,?) — ^(t/i)T^Therefore, on link 1, there is no wave; 
and on link 2, there is a forward shock wave when Dj < D2, a forward rarefaction 
wave when Di > D2, and no wave when Di = D2- 

Case 2 Link 1 is OC, and link 2 is UC. That is, 5i < Di == Ci and D2 < S2 = €2- In 
this case, Di — S2 — C2- From TheoremBJ we have that U{^ — U\ (0^ , t) = U2 = 
1/2(0^,0 — (Ci,C2) and ^1^2(0,?) — Ci. Therefore, on link 1, there is a back- 
ward rarefaction wave when Si < C\ and no wave when S\ =C\; and on link 2, 
there is a forward rarefaction wave when D2 < S2 and no wave when D2 =82- 

Case 3 Link 1 is OC, and link 2 is SOC. That is. Si < Di = C\ and ^2 < £>2 = C2. In 
this case, Ci =Di > S2- From TheoremB we have that U{^ — U\ (0^ , t) — U2 = 
U2{0^,t) — U2 and ^1^2 (0,f) = q{U2)- Therefore, on link 1 there is a backward 
shock wave when S\ > S2, a backward rarefaction wave when Si < S2, and no 
wave when Si = S2\ and on link 2, there is no wave. 

Case 4 Link 1 is SUC, and link 2 is OC, and q{Ui) < q{U2)- That is, £>i < ^2 < £>2 = 
Si =Ci. From Theorem[3]we have that f/f ^Ui{Q-,t) ^U^ =t/2(0+,f) = Ui 
and (71^2(0,0 = q{Ui). Therefore, on link 1, there is no wave; and on link 2, 
there is a forward shock wave. 

Case 5 Link 1 is SUC, and link 2 is SOC, and q{Ui) > ^(t/2). That is, ^2 < £>i < 5i = 
D2=Ci. From Theoremjsjwe have that f/f = t/i(0",f) = U^ = U2{0+,t) = U2 
and ^1^2(0,0 = q{U2)- Therefore, on link 1, there is a backward shock wave; 
and on link 2, there is a no wave. 

Case 6 Link 1 is SUC, and link 2 is SOC, and q{Ui) = ^(f/2). That is, Di = ^2 < £>2 == 
^i = Ci. From Theoremplwe have that C/f — U2 — iDi,S2), and ^1^2(0,?) = 
^(L'^i) — q{U2)- Therefore, on link 1, there is no wave; and on link 2, there 
is no wave. In this case, there can exist interior states on link 1 or link 2: 
Ui{Q-,t) = t/f and min{D2(0+, 0,52(0+,?)} > £>! or 1/2(0+, f) = U:^ and 
min{£)i(0-,r),5i(0-,0}>£>i. 

Obviously the results of stationary states and kinematic waves above are consistent 
with those in (| Lebacque|1996 ). That is, the new solution framework yields the same 
wave solutions as traditional approaches. 

The solutions of stationary states for the six cases are also shown in Figure l7J 
where figures (a)-(f) are for cases 1-6 respectively. In these figures, both the up- 
stream and downstream links share the same supply-demand diagram. From ini- 
tial conditions Ui and U2 we can first draw the pair (Di,52), from which we can 
determine upstream and downstream stationary states accordingly. Further we can 
summarize the solutions of stationary states in Figure [slin the (Di,5'2) space. This 
figure also demonstrates a graphical scheme for solving the stationary states as fol- 
lows. First, from initial Ui we draw a vertical line (thin pink line with an arrow), 
from initial U2 we draw a horizontal line (thin pink line with an arrow), and the 
intersection point is (Di,5'2). Then, if the intersection point is above the line OA, 
we draw a vertical line (thick blue line with arrow), and its intersection with AC 
gives the stationary states; if the intersection point is below the line OA, we draw a 
horizontal line (thick blue line with arrow), and its intersection with AC gives the 
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Stationary states; if the intersection point is on the Une OA, we draw both a vertical 
line (thick blue line with arrow) and a horizontal line (thick blue line with arrow), 
and their intersections with AC are the stationary states for the upstream and the 
downstream links respectively. Note that this scheme also works when the upstream 
and downstream links do not have the same fundamental diagrams but the same 
supply-demand diagram, i.e., Ci = €2- 




U[ = u{ = t/j 



Fig. 8 Solution of stationary states for the Riemann problem for the LWR model 



5.2 The inhomogeneous LWR model 



When Ci 7^ C2, then the road is inhomogeneous, and there is a discontinuity in the 
fundamental diagram at x = 0. In Figure l9j we demonstrate a graphical scheme for 
solving stationary states in the (01,82) space for the inhomogeneous LWR model. 
We take Figure |9ja) as an example, in which Ci < C2. First, from initial Ui we draw 
a vertical line (thin pink line with an arrow), from initial U2 we draw a horizontal 
line (thin pink line with an arrow), and the intersection point is (Di,5'2). Then, if 
the intersection point is above the line OA, we draw a vertical line (thick blue line 
with arrow), and its intersections with ACi and AC2 are the stationary states on 
links 1 and 2 respectively; if the intersection point (01,82) is below the line OA, we 
draw a horizontal line (thick blue line with arrow), and its intersections with ACi 
and AC2 are the stationary states on links 1 and 2 respectively; if the intersection 
point (01,82) is on the Une OA, we draw both a vertical line (thick blue Une with 
arrow) and a horizontal line (thick blue line with arrow), and their intersections with 
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ACi and AC2 are the stationary states for the upstream and the downstream links 
respectively. This scheme is the same as that for the homogeneous LWR model. 
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Fig. 9 Solution of stationary states of the Riemann problem for the inhomogeneous LWR model 



In ( Jin and Zhang|[2C)03a| l, the Riemann problem for the inhomogeneous LWR 
model was solved as a resonant nonlinear system, and ten types of wave solu- 
tions were obtained. For example, wave solutions of Type 1 can be obtained in the 
new framework as follows. Both Ui and U2 are UC, D2 < Di <C2 — S2, and Ci 
may be greater or smaller than €2- From Figure [9] or Theorem [3] we can see that 
U[' — (Di,Ci) — Ui, U2 ~ {Di,C2), there is no wave on link 1, there is a forward 
rarefaction wave on link 2, and ^'1^2(0,?) = q{U\)- It is easy to check that the wave 
solutions of other types are also consistent. 



6 Asymptotic traffic dynamics on an inhomogeneous ring road 



In this section we consider the inhomogeneous ring road with length L shown in 
Figure 10 in which the traffic direction is shown by the arrow. The ring road is 
composed of two homogeneous links: link 1 with capacity C\ for x € [0,Li], link 
2 with capacity C2 for x e [L\,L], the upstream boundary of link 1 is denoted as 
boundary 1, and the downstream boundary as boundary 2. Here we assume that 
link 1 is a bottleneck; i.e., C\ < €2- For example, such a bottleneck can be caused 
by a smaller number of lanes. We assume the fundamental relationships for two 
links as ^ = Qi{p) and p =^1(7) for;c e [0,Li], and q = Q2{p) and p = R2{y) for 
x€ [Li,L]. 
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boundary'2 



Fig. 10 A ling road 

6.1 Asymptotic stationary and interior states 

When the ring road reaches asymptotic stationary states, the flux at any location 
is the same, e.g., q. As we know, the asymptotic stationary state on a link can be 
uniformly UC, uniformly SOC, or a stationary shock wave (SS) connecting an up- 
stream sue state and a downstream SOC state (Bultelle et al.|1998). Then all pos- 



sible combinations of stationary states are listed in Table [T] and explained in the 
following. 

• When link 1 is UC at {q,C\) with <? < Ci, we have the following scenarios, (a) 
From Theoremlslit is possible that link 2 is UC at (g',C2), and the total number 
of vehicles on the ring road is A^„ ^ R\ {q/C\ )L\ + R2 {q/C2){L — Li). (b) If link 
2 is SS with upstream {q,C2) and downstream {C2,q), we have that q = Ci and 
link 1 is critical at (Ci,Ci). Assuming that link 2 is SUC for x e [^-1,^2] and 
SOC for X e [i2,i']- In this case, the total number of vehicles on the ring road is 
Nb=Riil)Li+R2{Ci/C2)iL2-Li)+R2iC2/Ci){L-L2).(c)lflmk2isSOC 
at {C2,q), we have from Theorem [3] that q = S2 =Ci at boundary 1. That is, link 
1 is critical at (Ci , Ci ) . In this case, the total number of vehicles on the ring road 
isNc=Ri{l)Li+R2{C2/Ci){L-Li). 

• When link 1 is SOC at {Ci,q) with q <Ci, we have the following scenarios, (d) 
It is possible that link 2 is SOC at {C2,q), and the total number of vehicles on the 
ring road is N,i = Ri {C\/q)Li +R2{C2/q){L- U). If link 2 is UC at {q,C2), we 
have from Theorem l3] that q = C\ at boundary 2. If link 2 is SS with upstream 
{q,C2), we have from Theorem [3] that q ^ C\ at boundary 2. Thus these two 
scenarios are impossible, since q = C\ contradicts q <Ci. 

• When link 1 is SS with upstream {q,Ci) and downstream (Ci , q) with ^ < Ci , we 
have the following scenarios. If link 2 is UC at {q,C2), we have from Theoremls] 
that q = min{Ci ,C2} = C\ at boundary 2; if link 2 is SS with upstream {q,C2), 
we have from Theorem [3] that q = minjCi ,C2} = Ci at boundary 2; if link 2 is 



Supply-demand diagrams and the inhomogeneous LWR model 



23 



SOC at (C2 ,q),we have from Theoremlslthat q = min{C2 ,Ci} — Ci at boundary 
1 . All three of these scenarios are impossible, since q — Ci contradicts q <Ci. 



^^^ Link 2 
Link 1^^^.^ 


UC {q,C2) 


SS {q,C2) ^ (C2,q) 


SOC (C2,<?) 


UC(./,Ci) 


(a) 


(b) 


(c) 


SOC(C,,</) 


X 


X 


(d) 


SS{q.C)^{Cuq) 


X 


X 


X 



Table 1 All possible stationary states on the ring road in Figure 



& 



For all four scenarios of asymptotic stationary states on the ring road, different 
scenarios have a different number of vehicles since 

0<Na<Ri{l)Li+R2{Ci/C2){L-Li)<Nb<Ne^ 

n{l)Li+R2{C2/Ci){L-Li) <Nd <RiHLi+R2H{L-Li). 

Due to traffic conservation on the ring road, we can therefore determine the fi- 
nal stationary states by the initial number of vehicles A^ on the road as follows: 
(a) When N < Ri{l)Li +/?2(Ci/C2)(L-Li), links 1 and 2 will be asymptotically 
stationary at UC with {q,Ci) and {q,C2) respectively, where q is the solution of 
Ri{q/Ci)Li +R2iq/C2)iL~Li)^N; (b) When /?i(l)Li +7J2(Ci/C2)(L-Li) < 
A' < Ri{l)Li +R2{C2/Ci){L — Li), link 1 will be asymptotically stationary at criti- 
cal with (Ci,Ci), and link 2 at SS with (Ci,C2) forx e [^1,^2] and (C2,Ci) forx e 
(L2,i], where L2 is the solution of /?i (l)Li +/?2(Ci/C2)(L2 -ii) +/?2(C2/Ci)(L- 
L2) = A^; (c) When A^ = Ri{l)Li +R2{C2/Ci){L- Ly), link 1 will be asymptot- 
ically stationary at critical with (Ci,Ci), and link 2 at SOC with (C2,Ci); (d) 
When A^ > R\{l)Li +7?2(C2/Ci)(L - Li), links 1 and 2 will be asymptotically 
stationary at SOC with {C\,q) and {C2,q) respectively, where q is the solution of 

Rl{Ci/q)Lx +R2{C2/q){L-U) ^N. 

From Theorem [3] an interior state can occur at a boundary when its upstream 
demand equals the downstream supply, and its flux cannot be smaller than the de- 
mand or supply. In the following we consider possible asymptotic interior states on 



the ring road in Figure 10 First, at any location inside a uniform traffic stream on 
a homogeneous road, it is not possible to have interior states, since the upstream 
and downstream states are exactly the same at (Z),5) and D = S if and only if the 
traffic is D = 5 = C, in which case the interior states have to be the same as the 
stationary states. Thus, interior states can only exist around the interface between 
two uniform traffic streams when the upstream demand equals the downstream sup- 
ply, and we examine possible interior states in all four scenarios as follows, (a) The 
necessary condition for an interior state to exist at boundary 1 is (7 = Ci, i.e., when 
A^ = R\{\)Li +R2{Ci/C2){L — Li). From Theorem l3] the interior state can only ex- 
ist at X ~0^, but not x = 0+. The necessary condition for an interior state to exist 
at boundary 2 is ^ = C2, which is not possible, (b) It is not possible for interior 
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States to exist at either boundary 1 or 2, but it is possible for an interior state to exist 
around the SS interface at x = L2- From Theorem l3] the interior state can only exist 
at X = L2 OT X = L2- (c) It is not possible for an interior state to exist at boundary 
1, but it is possible for an interior state to exist around boundary 2. That is, when 
N = Ri{l )Li +^2 (Ci /C2) (i — Li), from Theoremis] the interior state can only exist 
at X — L^, but not x — L^ . (d) It is not possible for interior states to exist at either 
boundary 1 or 2. 

In summary, there can exist three types of interior states: (a) When N = Ni = 
Ri{l)Li + R2{Ci /C2){L — Li), an interior state can only exist at x := 0^ = L^; 
(b) When 7;i(l)Li +7;2(Ci/C2)(L-Li) < A^ <7;i(l)Li +7;2(C2/Ci)(L-Li), an 
interior state can exist at x = L2 or x = LJ; (c) When N = N^ = Ri{l)Li + 
R2{C2/Ci){L — Li), an interior state can only exist atx — L^. 



6.2 Numerical examples 

In this subsection, we study asymptotic traffic dynamics on the inhomogeneous ring 



road in Figure 10 with L = 6001 = 16.8 km, Lj = 100/=2.8 km, and the location 
dependent speed-density relationships are based on ( Kerner and Konhauser||I994 
[Herrmann and Kerner|1998l ) 



y(p,a(x)) =5.0461 



l+exp{[— ^-^ 0.25]/0.06} ) -3.72x10 



-1 

6 



l/T, 



where the relaxation time T = 5 s; the unit length / = 0.028 km; the free flow speed 
Vf = 27.8 m/s; the jam density of a single lane Pj = 180 veh/km/lane. Here the 
number of lanes a{x) = 1 for link 1 and a{x) — 2 for link 2. The corresponding 
fundamental diagram q = Q{p,a{x)) is non-convex but unimodal in density p. In 
addition, Ci = 0.7091 veh/s, and C2 = 2Ci. Thus we can compute /?i(l)=35.8944 
veh/km, RiiCi/Ci) = /?2(2)=26.4162 veh/km, and RiiCi/Q) = /;2(2)=1 18.3550 
veh/km. Hence Ni = /?i(l)Li +/?2(Ci/C2)(L - Li) = 470.3311 veh, and A^3 = 
Ri{l)Li +R2{C2/Ci){L-Li) = 1757.4746 veh. 
Here we consider the following initial condition: 

p(x,0)=fl(x)(po + 3sin?P),xe[0,L], 

v(x,0) =y(p(x,0),fl(x)), xe[0,L]. ^ ^ 

Then, the total number of vehicles on the ring road is 

/■lOO/ 2nx 450 

N = 2poL- / (po + 3sin )dx =11001 po 1 

Jo L K 

When Po = 15.4007 veh/km, N = N\ and we observe an interior state at x = 0^ = 
L ; when po = 57.1911 veh/km, N = Nj, and we observe an interior state at x = Lj^ ; 
when Po G (15.4007,57.1911), we observe an interior state at x = Lj or x = LJ, 
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where L2 is the solution of /?i(l)Li +R2{\)[L2- Li) +R2{2){L- L2) ^N. For 
example, when po = 28 veh/km, A^ ~ 858.3893 veh, and 

■1^ 



L2 = 



N-iRi{l)-R2ii))Li-R2i2)L 



= 449.2561/. 



0.8 

d- "-^ 

p 0.4 

=^ 0.2 





100 200 300 400 500 600 



100 200 
^(100.0625, 1.3722) 



(449.1875,0.9684 



300 
x/l 



(599.9375, 1.0656). 



Fig. 11 Solutions of p, v, and q atT = 24000 s for initial conditions in |3l| : solid lines with stars 
for Po = 15.4007 veh/km, dashed lines with circles for po = 28 veh/km, and dash-dotted lines for 
po = 57.1911. 



In the following, we simulate traffic dynamics on the ring road for three different 
initial po: 15.4007 veh/km, 28 veh/km, and 57.1911 veh/km. Here we use the Go- 



dunov finite difference equation in ( 18 i and the supply-demand method in ( 19 1 for 
computing boundary fluxes. The simulation time is T = 4800T = 24000 s. We parti- 
tion the road [0,L] into A^ = 4800 cells and the time interval [0, T] into K = 240000 
steps. Hence, the length of each cell is Ax = 3.5 m and the leng th of each time step 



is At = 0.1 s. The CFL condition number (Courant et al. 



1928[)is y/^ <0.79<1. 
where the bottom 



The results for the three initial conditions are shown in Figure 
figure shows the locations and fluxes of all three interior states. From the figure, we 
can see that each of the three interior states only exists in one cell, and the locations 
of interior states are exactly as predicted above. Note that the top right figure of 



Figure 18 in (Jin and Zhang 2003a) also demonstrates the existence of an interior 
state, which is at the interface of a stationary shock. 
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7 Conclusion 

In this paper we first reviewed the definitions of the supply and demand functions 
and the discrete supply-demand method for computing boundary fluxes. We then in- 
troduced the supply-demand diagram of a roadway and a new framework for solving 
the Riemann problem of the inhomogeneous LWR model in supply-demand space. 
In this framework, each fink can have asymptotic interior and stationary states near 
the boundary, and the wave on each link is determined by the Riemann problem 
of the homogeneous LWR model with stationary and initial states for initial condi- 
tions. We have derived conditions for admissible stationary and interior states and 
introduced an entropy condition based on the discrete supply-demand method for 
computing boundary fluxes. We then proved that solutions to the Riemann problem 
exist and are unique and demonstrated that these solutions are consistent with those 
in literature for both the homogeneous and inhomogeneous LWR models. We also 
presented a graphical approach for finding the asymptotic stationary states with the 
help of supply-demand diagrams. Finally, we discussed the asymptotic stationary 
states on a ring road with arbitrary initial conditions and demonstrated with numer- 
ical examples that the existence and properties of the interior states are as predicted 
in this framework. 

Unlike existing studies of the homogeneous or inhomogeneous LWR models, 
this study analyzes traffic dynamics in supply-demand space. In this framework, the 
discrete supply-demand method is applied as an entropy condition. In this sense, our 
study provides a new approach for constructing convergent solutions of finite dif- 



ference equations arising in a Godunov method ( 18 i with a supply-demand method 
( [T9| for both the homogeneous and inhomogeneous LWR models. We have demon- 
strated that this new approach can successfully predict the existence and properties 
of interior states in numerical solutions. However, note that interior states take only 
one cell in numerical solutions and vanish as we diminish the cell size. In this sense, 
the interior states are inconsequential to solutions of the Riemann problem. 

Compared with existing studies, the new approach in the supply-demand frame- 
work is much simpler. In addition, since supply-demand methods have been pro- 



posed for computing fluxes through other junctions in general road networks (Da- 



ganzo 1995a Lebacque 1996 Jin and Zhang 2003b; Jin 2003 1, our framework could 



be extended to constructing solutions to the Riemann problem in these models. In 
(Jin||2009]l, we successfully applied this framework to analyze the Riemann prob- 



lem of merging traffic flow. In addition, one could also apply this new framework to 
analyze asymptotic traffic dynamics in a road network, such as the diverge-merge 
network studied in pin|2008| ). 
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